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FREFACE 

The  results  reported  here  are  derived  from  a  continuing 
theoretical  study  of  light  scattering  in  atmospheres,  supported  by 
the  U-S.  Air  Force  under  Project  RAND.  The  results  have  application, 
not  only  in  the  investigation  of  atmospheres  of  the  earth  and  other 
planets,  but  also  in  colloid  research.  This  publication  has  been 
prepared  for  presentation  at  the  Interdisciplinary  Conference  on 
Electromagnetic  Scattering,  Potsdam,  New  York,  August  19^2,  sponsored 
by  The  American  Chemical  Society,  and  by  the  (J.S.  Air  Force  Cambridge 


Research  Laboratories. 
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SUMMARY 


The  angular  scattering  and  polarization  properties  of 
polydispersed  suspensions  of  nonabsorbing  and  partially  absorbing 
spheres  have  been  computed  using  the  complete  Mie  series.  The 
results  with  three  types  of  size  distributions  are  presented  and 
compared  with  observations.  These  show  a  strong  dependence  of 
angular  intensity  and  polarization  patterns  on  the  size  distribu¬ 
tion,  the  size  range,  and  the  dielectric  and  absorbing  properties 
of  the  individual  particles.  A  peculiarity  of  scattering  at  angles 
near  45°,  observed  experimentally  and  independently  by  two  authors, 
is  corroborated  by  the  numerical  results.  Prominent  observational 
features  characteristic  of  natural  fog,  such  as  an  extremely  bright 
aureole,  rainbows  and  counter  coronas  are  reproduced  in  a  model 
corresponding  to  a  cloud  of  spherical  water  droplets,  with  a  wide 
distribution  in  droplet  radius  and  a  maximum  concentration  at  a 
4-micron  radius. 
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j 

Scattering  and  Polarization  Properties  of  Polydispersed 

* 

Suspensions  with  Partial  Absorption 


Introduction 

In  studying  the  theoretical  scattering  properties  of  poly¬ 
dispersed  suspensions,  such  as  natural  hazes  and  clouds  in  the 
atmosphere,  two  important  facts  should  be  considered:  (a)  Because 
of  the  finite  size  of  the  particles,  even  if  these  are  assumed 
spherical,  the  angular  scattering  characteristics  cannot  be 
approximated  with  sufficient  accuracy  by  asymptotic  expressions 
based  on  geometric  optics  or  Green's  function  approximations  for 
the  internal  field.  The  complete  Mie  series  must  be  used  for  each 
particle,  (b)  The  size  distribution  function  used  must  reproduce 
as  faithfully  as  possible  the  size  spectrum  of  the  suspension  in 
question. 

In  this  paper  we  shall  present  briefly  examples  simulating 
atmospheric  water  clouds  and  hazes  with  size  distributions  that  are 
sufficiently  realistic  to  allow  comparisons  with  observational  data. 
The  capabilities  of  modern  high-speed  computers  in  generating  single¬ 
particle  Mie  scattering  functions,  both  with  real  and  complex  index 
of  refraction,  has  been  adequately  described  elsewhere  (Deirmendjian, 
et  al,  1961;  Giese,  1961) .  Therefore  we  shall  describe  only  the 

* 

Parts  of  this  paper  were  presented  at  the  International  Symposium 
on  Radiation  of  the  Radiation  Commission,  I.A.M.A.P.,  Vienna, 

August  1961. 
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size-distributicn  functions  used  and  Justify  their  choice. 

Size -distribution  function 

The  function  chosen  is  a  generalization  of  that  first  proposed 
by  Khrgian  and  Mazin  (1952,  195^)  for  clouds.  A  survey  of  various 
proposed  functions  showed  that  the  latter  type  seems  to  fit  fairly 
well  with  measurements  both  on  natural  water  clouds  and  aerosols,  and 

e 

it  has  the  great  advantage  that  its  parameters  have  readily  interpretable 
physical  meaning.  In  its  most  general  form  this  function  may  be 
written 

/  .  OL  -br^  /  _ . 

n(r)  =  a  r  e  (1) 

where  n(r)  is  the  volume  concentration  at  the  radius  r  and  a,  a,  b  and 
7  are  positive  constants.  For  a  particular  choice  of  a  and  7,  the 
remaining  constants  can  be  uniquely  determined  by  the  total  number  of 
particles  N  per  unit  volume,  and  by  the  critical  or  mode  radius  rc, 
where  the  concentration  is  at  a  maximum.  Thus  the  constant  a  can  be 
determined  from  the  integral 

,00  OH-1 

N  =  \  n(r)dr  *  2.  b  V  (~i)  ,  (2) 

and  b  from  the  derivative,  which  can  be  written 

y 

—  n(r)  =  a  r  e  (a  -  7br' )  (3) 

which  vanishes  when 
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b  = 


a 


—  ,  r  =  r 

7  c 


7  r 

as  well  as  for  r  =  0  and  oo  . 

In  particular,  with  the  choice  CH 
100  cm  as  a  nominal  concentration, 
function 


=  6,  7  =  1, 
one  gets  the 


(4) 


r__  =  Up ,  and  N  = 
di stribution 


n(r) 


2.373  r6  e~1-5r 


cm 


•3 


(5) 


when  the  radius  is  expressed  in  microns  and  the  concentration  per 
cubic  centimeter.  This  model,  marked  "cloud"  in  Fig.  1, 
reproduces  fairly  well  some  of  the  distributions  found  by  Durbin  (1959) 
for  cumulous  clouds  having  a  depth  between  230  and  2100  m,  as  averaged 
and  reproduced  as  "cumulus  Type  I"  in  a  diagram  by  Singleton  and 
Smith  (i960). 

The  discontinuous  distribution  called  "haze  C"  in  Fig.  1  roughly 
approximates  the  continental  hazes  measured  by  various  workers  and 
used  by  the  author  previously  (Deirmendjian  1957,  1959,  i960). 

The  other  continuous  distribution  in  Fig.  1,  marked  "haze  M"  is 
obtained  by  putting  CC  =  1,  7  =  l/2  in  (l)  with  r^  =  resulting 

in  the  specific  function 


n(r) 


5-33  x  10 


4 


re 


8.9I+4  \f r 


(6) 


This  was  found  to  agree  with  unpublished  counts  of  aerosol  particles 
made  in  Los  Angeles  by  J .  Y.  Gilbert  (1954),  and  may  be  taken  to 
represent  coastal  conditions.  The  differences  from  the  model  C  haze_ 
are  considerable,  in  that  there  is  a  much  larger  concentration  of 


Fig. I  —  Three  size-distribution  functions  used  in  the 


integration  of  the  Mie  functions.  Total  concentration  IOO  cm 
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particles  in  the  0.15- to  3 .5-micron  range.  The  distribution  (6)  seems 
also  to  correspond  fairly  well  to  that  of  stratospheric  aerosols 
found  by  Junge,  Chagnon  and  Manson  (l9<5l),  at  least  in  the  range 
0 . 2  <  r  <l(i . 

The  logarithmic  representation  in  Fig.  1  is  chosen  for  convenience 
only  and  does  not  give  a  good  idea  of  the  distribution.  A  linear 
representation  of  the  continuous  distributions  would  result  in  bell¬ 
shaped  curves  which  are  skewed  toward  the  larger  sizes.  However, 
since  the  integrals  under  all  three  distributions  in  Fig.  1  give 
100  particles  per  cubic  cm,  the  curves  do  provide  an  idea  of  the 
relative  amounts  in  each  size  interval. 

Integration  of  the  Mie  functions  with  size  distribution 

To  get  the  total  scattering  and  absorption  cross-sections  per 
unit  volume  of  the  cloud  or  haze, one  has  to  integrate  the  Mie  single¬ 
particle  cross-sections  with  respect  to  the  appropriate  distribution. 
Making  the  convenient  change  of  variable  r  =  k  "'"x,  where  x  is  the 
Mie  size  parameter  and  k  =  2rt/a  is  the  free-space  propagation 
constant  at  the  wavelength  J\,  one  can  rewrite  the  size  distribution 
function  (l)  in  the  form 


f(x,k)  =  a  k”^"a  x“  exp  ^  -  b  (x/k)^ 


(7) 


-3 

where  the  factor  k  is  introduced  for  convenience.  It  is  easy  to 
show  then  that  the  volume  cross-sections  p  may  be  obtained  by  the 
integration 


p(m,^;  x1,x2)  =  « 


2 


x 


x 


1 


f(x,k)  K(m,x)  dx 


(8) 
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where  m  is  the  complex  or  the  real  index  of  refraction  of  the  particles, 
K  is  the  Mie  extinction,  scattering,  or  absorption  cross-section,  as 
the  case  may  be,  divided  by  the  geometrical  cross-section  (Deirmendj ian 
et  al.,  1961),  and  x^  and  Xg  correspond  to  the  lower  and  upper  limits  in 
the  particle  size. 

The  elements  of  the  normalized  scattering  matrix  or  phase  matrix 
as  defined  by  Chandrasekhar  (1950),  are  given  by  the  integrals 


^  (m,6,y;x1,x2)  =  i  \  f(x,k)  ij(m,x,e)dx,  J  =  1,2, 3, k  (9) 


where  9  is  the  scattering  angle,  and  i.(m,x,9)  are  the  dimensionless 
Mie  intensity  functions  computed  previously  (Deirmendjian,  1961), 
related  to  the  intensity  or  phase  matrix  operating  on  the  Stokes  vector 
of  the  incident  radiation  in  the  manner  indicated  by  various  authors 
(Sekera,  1957;  van  de  Hulst,  1957) • 

It  is  in  the  behavior  of  the  integrals  (8)  and  (9)  as  a  function 
of  Xg  that  the  advantage  of  the  distribution  functions  of  the  type  (l) 
becomes  evident.  Numerical  tests  using  our  IBM7090  program  show,  in 
fact,  that  in  this  case,  these  integrals  always  converge  to  a  constant 
value  for  all  x0^>  X,  where  X  is  a  finite  number  depending  on  the 
logarithmic  derivative  of  the  function  n(r)  with  respect  to  r.  In 
particular,  it  appears  that  a  sufficient  condition  for  convergence 
to  begin,  for  all  integrals  (8),  (9)  including  the  forward  scattering 
case  with  9=0,  is  that  (or  correspondingly  Xg)  be  such  that 


Toi"T  log  n<r>  =  a  (r)  -  1  > 
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Solving  the  inequality  to  the  right,  we  get  the  condition  for 
convergence  to  begin 

r2>rc  t1  +  h/a)^7  (11) 

This  is  a  very  convenient  property  in  regard  to  the  contribution  of 
the  larger  particles  in  a  suspension,  where  the  corresponding  Mie 
series  are  very  long  and  machine  computation  time  becomes  costly. 

It  is  also  a  physically  justifiable  feature  of  long-lived  atmospheric 
aerosol  and  cloud  particles,  which  are  naturally  limited  in  the 
maximum  size.  Distributions  such  as  model  C  with  n(r)~  r  ^  are  not 
convenient  because  the  integrated  volume  extinction^  and  scattering 
in  the  forward  area  depend  on  the  choice  of  the  upper  limit  r^. 

The  numerical  results  further  show  that,  for  x^  in  the  region 
defined  by  (ll),  the  integrals  (8)  and  (9),  when  plotted  as  a 
function  of  x^,  are  smooth  and  monotonic  curves,  with  little  change 
in  value  when  the  integration  interval  Ax  is  varied.  This  is  a 
remarkable  property,  connected  with  the  choice  of  continuous  size 
distribution  functions  of  the  type  (l),  since  it  is  known 
(Penndorf,  1962;  Deirmendjian,  et  al,  1961) ,  that  the  intensity  for 
single  particles,  especially  in  the  backward  hemisphere,  shows  large 
fluctuations  as  a  function  of  the  size,  whose  amplitude  and 
frequency  increase  with  x.  It  means  that  the  angular  scattering 
pattern  of  polydispersed  suspensions,  with  a  more  or  less  wide  and 
continuous  size  spectrum,  should  be  a  smooth  function  of  0,  without 
the  numerous  coronas,  rainbows  and  countercoronas  characteristic 
of  single  particles  or  monodispersed  suspensions. 


It  should  be  noted  from  the  form  of  the  integral  (9)  that  the 
phase  functions  are  independent  of  the  total  concentration  N  of 
particles  per  unit  volume,  as  should  be  evident  from  physical 
considerations.  This  statement  is  of  course  valid  only  when  the 
concentration  is  small  enough  that  only  incoherent  or  independent 
scattering  need  be  considered. 

Volume  extinction  and  albedo 

Examples  of  the  extinction  per  unit  path  for  the  three  models 
and  at  selected  wavelengths  are  shown  in  Table  1.  The  particles 
are  assumed  to  be  made  of  liquid  water,  with  the  real  or  the  complex 
index  of  refraction  indicated  in  Table  1,  adapted  from  Centeno  (19^1) 
The  table  also  lists  the  computed  albedo  A  per  unit  volume  of  space, 
that  is,  the  fraction  of  energy  scattered  in  all  directions  weighted 
by  the  distribution  function,  compared  to  that  scattered  and  absorbed 
The  numbers  in  Table  1  illustrate  clearly  the  effects  of  changes  in 
distribution  function,  in  the  relative  size,  and  in  the  index  of 
refraction  of  the  particles.  Noted  particularly  are  the  similarities 
and  differences  between  the  two  haze  models.  The  lower  albedo  for 
model  C  at  93.07m.  and  -96.05m  is  explained  by  the  larger  proportion 
of  small  particles  in  it  as  compared  to  model  M  (see  Fig.  l),  such 
particles  being  relatively  more  efficient  real  absorbers  than  the 
larger  ones  for  any  given  complex  index.  These  differences  are 
even  more  pronounced  when  comparing  the  haze  and  cloud  distributions, 
especially  at  JOLOm  • 

The  extinction  values  for  the  cloud  model  show  the  interesting 
fact  that  these  should  increase  with  wavelength  in  the  near-infrared 
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region,  even  though  the  relative  size  decreases,  partly  because  of 

the  increase  in  the  absorption  coefficient  of  liquid  water. 

A  comparison  of  the  corresponding  extinction  coefficients 

under  Table  1  for  the  model  C  haze,  with  those  approximated 

previously  by  Deirmendjian  (i960)  by  means  of  an  empirical  extension 

of  van  de  Hulst's  (1957)  analytical  approximation,  shows  very  good 

agreement,  indicating  the  sound  basis  of  van  de  Hulst's  original 

analysis.  Those  corresponding  to  a  cloud  also  agree  in  order  of 

magnitude  and  wavelength  dependence,  but  not  in  detail,  because 

a  different  model,  with  r  =  3-5u  ant  a  narrower  symmetrical 

c 

distribution,  was  used  in  the  earlier  approximate  calculation  by 
Deirmendj  ian  ( i960 ) . 

Normalized  intensity  functions  P.^  and  Pg 

Figures  2  to  9  represent  examples  of  the  integrated  intensity 
functions  of  phase  matrix  elements  (solid  line)  and  Pg/Ujt 

(dashed  lines)  which  are  normalized  so  that  they  satisfy  the 
condition: 

^  +  p2(0^  dco  =  1  (12) 

A 

The  ordinates  in  all  these  figures  are  on  a  logarithmic  scale,  while 
the  abscissa  is  linear  in  the  scattering  angle  0  measured  from  the 
forward  direction.  For  unpolarized  and  parallel  monochromatic 
incident  energy  of  unit  flux,  the  degree  of  partial  linear 
polarization  is  given  by  the  ratio 

P1(0)  - 

P1(9)  +  P2(0) 
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with  positive  values  corresponding  to  maximum  electric-vector  amplitude 
perpendicular  to  the  scattering  plane.  The  maximum  numerical  value 
and  position  of  this  parameter  is  indicated  on  all  eight  figures. 

Figure  2  corresponds  to  a  model  C  water  haze  (maximum  particle 
radius  of  5ji)  illuminated  by  infared  radiation  at  ^5*3^*  The 
scattering  pattern  is  smooth  with  strong  anisotropy  equivalent  to  a 

p 

ratio  of  10  between  forward  and  backward  scattering.  A  maximum 
positive  polarization  of  0.71  occurs  at  0  =  100°,  explained  by  the 
small  relative  size  of  the  particles  and  their  absorption  properties. 

Figure  3  shows  the  scattering  pattern  for  the  same  haze  when 
illuminated  by  visible  light  at  0iO.45g,  where  the  refractive  index 
of  water  is  assumed  real  and  equal  to  1.34.  If  Fig.  3  is  compared 
to  Fig.  2,  it  is  seen  that  the  anisotropy  has  increased  by  one  order 
of  magnitude  and  the  maximum  polarization,  though  still  at  0  =  100°, 
has  decreased  in  value  to  0.49.  Both  changes  are  due  to  the  larger 
relative  size  of  the  particles.  The  pattern  for  J^0.70|i  radiation 
(not  shown)  is  similar  except  that  the  maximum  at  0  =  0  is  reduced 
by  about  62  per  cent  . 

Figures  4  and  5  correspond  to  the  situation  for  a  water  haze 
with  a  distribution  according  to  model  M,  at  the  wavelengths  of 
0.70  and  0  .b’^x  respectively.  The  intense  aureole  in  the  region 
0O<e<10°,  is  more  pronounced  at  the  shorter  wavelength  here  than 
in  model  C,  with  a  reduction  of  50  per  cent  atCX0.70n-  Another 
difference  is  in  the  polarization  pattern,  which  has  migrated 
toward  larger  angles  with  a  maximum  of  about  0.40  at  9  =  152  •  In 
the  region  20°<  0<1OO°,  the  polarization  is  weak  and  negative.  All 


HAZE  MODEL  C 

X5.30  fj.  A  =  0.80 


120  140  160  180 


9  ( DEGREES) 


Fig.  2 

Integrated  and  normalized  intensity  functions  P^/4-Jt  (solid  line) 
and  P  /4n  (dashed  line)  for  haze  particles  with  complex  index 
(see  Table  l)  and  infrared  illumination.  Computed  values  at 
0:  0(5)20(l0)l50(5)l80o  joined  by  straight  lines. 
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B (DEGREES) 


Fig.  3 

Same  as  Fig.  2  but  for  particles  with  real  index  and  visible 
light.  Computed  values  at  0:  0(l)5>10(l0)l50(5)l70(2)l80° . 
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Fig. 4 — Intensity  functions  for  haze  particles  at  \0.70/z  and  with 
real  index  1.33  but  a  different  distribution  than  in  Fig. 3 
Computed  values  at  9-  0(2.5)  20(10)  130(2.5)  180° 
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these  changes  in  going  from  a  model  C  to  a  model  M  distribution  are 
connected  with  the  relatively  larger  amount  of  particles  in  the 
range  0'15<Cr<£3'5\i  in  the  second  model  as  compared  to  the  first 
(see  Fig.  l) . 

In  general,  the  unit-volume  scattering  characteristics  by  the 
haze  particles  must  be  closely  related  to  the  deviations  of  the 
observed  degree  of  polarization  and  the  position  of  the  neutral 
points  in  the  sunlit  sky  from  those  predicted  by  Rayleigh  multiple 
scattering,  as  discussed  by  Sekera  (1957)  and  others.  More 
particularly,  it  can  be  shown  (DeirmendJ ian,  1959)  that  the  most 
prominent  deviations  should  appear  in  the  brightness  of  the  region 
around  the  sun  called  the  aureole,  which  can  be  easily  explained  by 
the  effect  of  primary  scattering  of  sunlight  on  haze  particles.  In 
fact,  the  normalized  intensity  function  PA(0)  for  a  mixture  of  air  and 
particles  may  be  written  in  the  fora 


PA«>) 


%  +  PM  PM(e> 


BR+  s« 


(13) 


where  the  subscripts  R  and  M  stand  for  Rayleigh  and  Mie  particles 
respectively.  Usually,  even  on  clear  days,  the  attenuation 
coefficients  PR  and  in  the  visible  range,  are  of  the  same  order 
of  magnitude,  whereas  for  small  scattering  angles  say  0<C5°> 
pM  —  100  PR,  where  PM  =  !4  (p1+p2)  and  PR  =  %  (l+cos20)  .  Thus,  at 
such  angles,  the  large  particle  effect  will  dominate  in  the 
brightness  (and  hence  also  in  the  polarization). 

In  this  connection,  it  is  important  to  introduce  the  correct 
gradient  and  intensity  of  the  aureole  in  estimates  of  the  volume 


scattering  coefficient  based  on  an  integration  of  the  angular  intensity 
over  all  solid  angles.  Extrapolations  of  the  angular  intensity  in 
the  region  0°  =  0<^1O°,  when  the  corresponding  data  are  unavailable, 
may  lead  to  an  incorrect  estimate  of  both  the  volume  scattering 
coefficient  and  the  anisotropy  between  forward  and  back  scattering, 
and  hence  also  of  the  wavelength  dependence  of  these  parameters 
(Gibbons  et  al.,  1961) . 

The  exact  integrations  of  the  Mie  functions  can  now  be  used  to 
check  a  previous  approximation  by  Deirmendj ian  (1957*  1959)  based  on  the 
W.K.B.  method.  In  that  paper,  the  Mie  phase  functions  were  under¬ 
estimated  by  a  factor  4rt  (1959>  equations  (15)  >  (l6)  and  Table  3) 
because  of  an  oversight  in  the  definition  of  the  amplitudes  of  the 
scattered  field.  Taking  this  into  account,  the  following  table 
gives  a  comparison  of  the  approximate  and  exact  calculations  for  a 
model  C  haze  at 


Phase 

function 

Px(e)> 

model  C 

haze 

10° 

0 

0?25 

i° 

3° 

6° 

W.K.B.  method 

562 

157 

59.3 

31-5 

18.3 

(Deirmendjian,  1959) 

Exact  Mie  functions 

74.7 

67.1 

38.1 

20.0 

14.7 

(this  paper) 

The  large  deviations  appearing  above  are  explained  by  the  incorrect 
integration  with  respect  to  a  model  C  distribution  rather  than  by  the 
W.K.B.  approximation  itself,  which  is  very  good  for  individual 
particles  at  small  scattering  angles  (Deirmendjian,  1957)*  In  fact, 
if  the  original  integration  of  the  diffracted  intensity  were  carried 
to  a  realistic  finite  upper  limit  of  the  particle  size,  instead  of  an 


-18- 


infinite  limit  as  in  the  Struve  integral  (Deirmendj ian,  1959,  p.  231 ), 
the  use  of  the  W.K.B.  approximation  should  result  in  very  good 
estimates  of  the  aureole  for  this  particular  model.  Furthermore,  if 
a  distribution  of  the  type  (l)  is  used  instead,  the  convergence 
property  mentioned  in  the  second  section  of  this  paper,  suggests  that 
in  the  diffraction  area,  definite  integrals  of  the  form 

QD 

2+a  -bx  r_  ,  .  2  . 

x  e  |^J^(xsinS)J  dx 

must  exist,  with  a  finite  value  even  for  9=0. 

Figures  6  to  9  present  f^ie  normalized  intensity  obtained  from  a 
cloud-drop  size  distribution  described  by  (5)-  Figure  6  corresponds 
to  infrared  illumination  at  JAlOp.  A  considerable  anisotropy  and 
strong  maximum  polarization  at  9  *  110°  are  evident,  even  for  this 
long  wavelength;  otherwise  the  angular  scattering  pattern  is  quite 
smooth.  The  next  case,  shown  in  Fig.  7,  can  be  compared  with  Fig.  2 
for  haze.  It  is  seen  that  in  the  cloud  model  a  maximum  positive 
polarization  occurs  at  9  =  152°,  in  a  broad,  secondary  maximum  in 
the  intensity,  corresponding  to  an  infrared  rainbow.  The  anisotropy 
is  high  with  a  ratio  of  about  10  between  forward  and  back  scatter. 

Going  towards  the  visible,  Figs.  8  and  9  show  the  very  interesting 
patterns  for  a  cloud  illuminated  by  1*0.70  and  !X0.45p  radiation, 
respectively.  Note  first  the  very  pronounced  aureole  or  diffraction 
peak  in  the  region  0°  =  9^-  5°,  which  is  twice  as  bright  in  the 
blue  as  in  the  red  near  0°.  This  undoubtedly  accounts  for  the  very 
bright  "silver  lining"  at  the  borders  of  certain  clouds  when  near  the 
sun.  Next,  in  the  region  20°<C  0<8o°,  the  intensities  fall 


sin  9 


Q  (DEGREES) 

Fig.  6 

Intensity  functions  for  a  water  cloud  illuminated  "by  MOp 
infrared  radiation.  Computed  values  at  9:  0(5)180  . 
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9  (DEGREES) 

Fig.  7  ® 


5ame  as  Fig.  6  but  at  A5- 3M- • 

values  at  9:  0(2.5) 30(l0)l4o(2.5)l80°. 


0  20  40  60  80  100  120  140  160 

9 (DEGREES) 

Fig.  8 

Same  as  Fig.  6  but  at  A0.70|i  and  real  index  1.33 
Computed  values  at  0:  0(l)l0(5)l30(2)l80  . 
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Fig. 9  —  Same  as  Fig. 6  but  at  X0.45/J.  and  real  index  1.34 
Computed  values  at  9  ■■  0(1)  15(5)  120(2)  132(1)  180° 
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exponentially  with  the  scattering  angle,  with  little  or  no  polarization 
change.  The  region  8o°<^  0<CllO°  corresponds  to  a  minimum  in  intensity 
and  change  in  the  sign  of  the  polarization. 

The  region  110 0  =  l8o°  contains  distinct  featuies  which 
deserve  discussion.  The  most  salient  of  these  is  of  course  the 
pronounced  and  relatively  narrow  maximum  around  0  =  143°  •  This  is 
the  fogbow  or  white  rainbow  mentioned  by  Minnaert  (1954,  p.  183) . 
Comparing  Figs.  8  and  9>  we  see  that  there  is  no  true  dispersion, 
the  blue  bow  maximum  corresponding  to  the  red  in  position.  However, 
the  blue  bow  is  narrower  and  more  intense  and  its  polarization 
exceeds  that  of  the  red  bow.  This  must  be  explained  in  terms  of  the 
larger  relative  size  of  the  cloud  drops  with  respect  to  the  wavelength. 
This  trend  indicates  that  the  polarization  in  the  true  rainbow  should 
be  even  greater.  Furthermore,  in  the  true  rainbow,  which  occurs  at 

about  0  =  138°,  the  condition  x  - ►  co  applies  to  all  wavelengths 

and  the  dispersion  must  be  mostly  due  to  the  ref rangibility  of 
water.  The  interesting  fact  that  emerges  from  the  present  study 
is  that  a  wide  size  distribution  in  cloud  droplets  does  give  rise  to 
a  single  principal  bow  so  that  all  the  secondaries  and  supernumeraries, 
which  appear  by  scattering  on  single  cloud  droplets,  are  virtually 
suppressed  by  superposition.  There  is  an  indication  of  a  weak 
secondary  bow  around  0  =  122°  at  ^0.45n  only. 

Similarly  the  appearance  of  a  counter  corona  of  radius  of  about 
2°  at^0.45n  and  4°  at^0.7CV  is  interesting  for  the  same  reason. 

The  author  has  actually  observed  this  visually  on  the  vertical  wall 
of  a  towering  cumulus  from  a  jet  airliner  at  about  30,000  feet,  with 
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the  sun  near  the  horizon.  Note  that  the  polarization  in  the  counter¬ 
corona  has  the  opposite  sign  from  that  of  the  rainbow.  No  direct 

coronas  were  found  in  this  model,  however,  because  of  the  wide  size 

© 

distribution  assumed. 

Comparisons  with  observations  and  other  work 

A  most  remarkable  corroboration  of  the  above  cloud  model  comes 
from  the  careful  angular  scattering  and  polarization  measurements 
made  at  night  by  Pritchard  and  Elliott  (i960)  with  their  recording 
polar  nephelometer  on  small  volumes  of  natural  fog  (their  Fig.  13) . 
These  measurements,  made  with  a  green  filter,  reproduce  quite 
faithfully  most  of  the  salient  features  appearing  in  Figs.  8  and 
including  the  intense  aureole,  the  exponential  decrease  and  other 
gradients,  position  of  the  minimum,  position  and  intensity  of  the 
rainbow  and  its  polarization  (the  components  W  and  HII  of  Pritchard 
and  Elliott,  i960,  are  proportional  to  the  elements  P^  and  P ^ 
respectively).  Apparently  it  was  not  possible  to  cover  the  range 
170  9  =  l8o°,  hence  the  counter  corona  and  backward  peak  are  not 

present .  Considering  that  the  observations  correspond  to  ground  fog, 
in  which  the  size  distribution  is  not  necessarily  that  of  a  cumulous 
cloud  as  in  our  model,  the  agreement  is  indeed  remarkable.  Also  the 
measurements  include  the  effects  of  molecular  scattering  which  must 
have  some  influence  on  the  polarization  around  9  =  100°  (see 
Equation  (13)  )  . 

The  measurements  on  clear  nights  by  Pritchard  and  Elliott  (i960) 
also  show  good  agreement  with  our  haze  models  except  in  the 
polarization,  whose  maximum  does  appear  at  9  =  100°  but  has  a  larger 
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value  than  shown  on  Figures  3  or  5>  for  example,  because  of  the  strong 
influence  of  molecular  scattering  at  this  angle. 

The  theoretical  curves  should,  also  be  compared  with  measurements 
of  the  angular  intensity  only,  in  the  range  l6°  =  0  4  164°  reported 
by  Barteneva  (i960),  whose  work,  however,  is  not  as  detailed  as  that 
of  Pritchard  and  Elliott. 

Integrations  somewhat  similar  to  ours,  but  not  normalized  with 
respect  to  the  integrated  scattering  cross-sections,  and  carried  out 
to  the  constant  upper  limit  of  x 2  ■  40,  with  an  interval  of  Ax  =  1, 
have  been  reported  by  Giese  (1961).  In  contrast,  for  example,  our 
integration  intervals  and  range  in  x  for  Figure  9  is 
0.25(0.25)6c(0.50)l60.  Furthermore  Giese  has  used  distribution 
functions  of  the  type  n(r)~r"a  with  a  =  2,  2-5,  3  appropriate  to 
zodiacal  particles .  It  is  interesting  to  note  that  samples  of  his 
results  (Giese,  1961,  Fig.  9&)  show  some  of  the  features  appearing 
in  our  cloud  model,  including  a  broad  rainbow  arcund  0  =  l4o°,  even 
though  his  size  range  corresponds  to  a  haze  model.  This  is  explained 
by  the  size  distribution  function  used  by  Giese,  which  is  considerably 
less  steep  than  in  our  models  C  and  M. 

One  very  interesting  observational  peculiarity  of  scattering 
on  polydisperse  Mie  particles  must  be  mentioned,  whose  theoretical 
corroboration  could  be  provided  only  by  computing  the  normalized 
intensities  as  in  the  present  case.  This  is  the  apparent  constancy 
of  the  quantity  (P1+P2)/8«  around  0  =  40°,  with  a  value  of  about 
0.1  regardless  of  the  size  distribution  function,  as  evidenced  by 
Figs.  4,  5,  8  and  9«  A  check  into  our  tabulated  values  shows  that 
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this  is  true  even  if  the  integration  is  stopped,  at  smaller  sizes, 
provided  these  exceed  the  first  resonance  peak  in  the  scattering 
cross-section.  This  behavior  suggests  that  for  small  samples  of 
a  polydisperse  suspension,  the  concentration  of  the  particles  is 
directly  proportional  to  the  measured  intensity  scattered  around 
9  ■  40° ,  regardless  of  the  size  distribution  and  maximum  size  of  the 
particles,  provided  the  distribution  is  continuous  and  of  the  type 
(1),  and  the  index  of  refraction  is  real  and  remains  unchanged. 

The  experimental  facts  are  provided  by  the  analyses  of  two  independent 
sets  of  observations  made  on  entirely  different  media,  to  wit,  those 
reported  by  Barteneva  and  Bashilov  (1961)  for  atmospheric  aerosols, 
and  those  of  Tyler  (1961)  for  hydrosols.  Both  papers  reveal  a  very 
good  linear  correlation  between  the  scattering  function  at  9  -  4-5° 
and  the  volume  scattering  coefficient  or  degree  of  turbidity.  This 
property  is  also  mentioned  in  a  very  interesting  review  paper  by 
G.  V.  Rosenberg  (i960). 

The  functions  P^  and 

The  functions  P^  and  P^  integrated  and  normalized  according  to 
(9)  have  to  do  with  a  rotation  of  the  plane  of  polarization  and 
introduction  of  elliptical  polarization  by  polydisperse  suspensions, 
when  illuminated  by  polarized  light.  These  will  not  be  discussed  in 
detail  here,  except  to  show  one  example  in  Fig.  10  corresponding  to 
Model  M  at3v0«70n.  The  functions  ¥^(9)/bn  (solid  line)  and  P^(®)/4« 
(dashed  line)  are  plotted  linearly  against  the  scattering  angle. 

For  comparison,  the  corresponding  Rayleigh  value  for  P ^B)/hn  is 
shown  by  the  dotted  line  (the  Rayleigh  value  for  P^  is  zero 
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Fig.lO —  Integrated  and  normalized  functions  (solid 

line)  and  F^/4  7r  (dashed  line)  for  haze  particles  at  \0.70/x. 

Comparable  values,  shown  by  solid  and  open  dots, 
respectively, for  a  different  size  distribution  and  refractive 
index,  adapted  from  Fraser  (1959,  p.  123, "model  D") 
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everywhere).  There  are  no  experimental  values  that  could  be  used  for 
comparison  here.  Comparable  values  computed  by  Fraser  (1959)  are 
shown  by  solid  and  open  dots  in  Fig.  10.  Although  the  latter  were 
computed  for  m  =  1.^0  and  a  different  size  distribution  and  range, 
there  is  general  agreement  in  the  order  of  magnitude,  the  gradient 
and  the  position  of  the  maxima.  Fraser’s  results  are  based  on  a  Legendre 
series  expansion  due  to  Sekera  (1957) >  using  the  Mie  scattering 
coefficients  a^  and  b  ,  which  are  independent  of  9.  The  general 
agreement  mentioned  above  is  proof  of  the  correctness  of  the  method. 
However,  the  possibility  of  accurately  generating  the  complex 
scattered  amplitudes  by  means  of  high-speed  computers  (Deirmendjian 
et  al,  1961),  allows  a  much  more  reliable  direct  integration  of 
these  parameters. 

Some  conclusions 

On  the  basis  of  this  limited  sample  of  our  results  and  comparison 
with  available  observations,  the  following  conclusions  may  be 
reached: 

a.  The  classical  electromagnetic  theory  of  scattering  on 
homogeneous  finite  spherical  particles,  applied  to  a  polydispersed 
suspension  with  a  continuous  size  distribution,  can  be  used  to 
reproduce  rather  accurately  the  observed  angular  intensity  and  linear 
polarization  of  samples  of  natural  water  clouds  and  hazes  illuminated 
by  visible  light. 

b.  No  single  particle  in  such  a  suspension  exists,  whose  angular 
scattering  and  polarization  properties  correspond  to  the  integrated 
effects  from  all  the  particles,  which  tend  to  suppress  by 
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superposition  the  extreme  angular  fluctuations  characteristic  of 
single  particles.  The  introduction  of  such  parameters  as  the  "mean 
radius,"  the  "mode  radius,"  the  "mean  volume  radius,"  etc.,  is  not 
helpful  in  determining  these  properties  for  aggregates  of  particles, 
except  perhaps  in  the  case  of  very  small  particles  relative  to  the 
wavelength.  Our  results  further  indicate  that  the  angular  volume 
scattering  properties  are  uniquely  determined  by  the  type  of  size 
distribution,  the  range  of  sizes,  and  the  dielectric  and  conducting 
properties  of  the  particles  in  the  suspension,  independently  of  the 
total  concentration.  Hence,  in  the  case  of  optically  thin  samples 
of  such  suspensions  (where  multiple  scattering  is  negligible),  these 
parameters  may  be  deduced  from  careful  measurements  of  the  scattering 
properties  in  as  many  directions  as  possible  and  at  various 
wavelengths . 

c.  The  existence  of  a  solution  to  the  analogous  general  problem 
in  the  case  of  optically  thick  media,  i.e.,  that  of  the  unique 
determination  of  the  type,  size,  size  distribution,  and  size  range 
of  the  particles  in  suspension,  from  a  knowledge  of  the  angular 
intensity  and  polarization  of  the  diffusely  transmitted  and  reflected 
radiation,  as  discussed  Dy  other  participants  in  this  symposium, 
remains  to  be  demonstrated.  However,  in  the  case  of  certain  suspensions 
such  as  water-droplet  clouds  in  the  atmosphere,  the  existence  of 
prominent  features  such  as  rainbows  and  counter  coronas  suggests  that 
some  information  on  the  size  and  size  distribution  of  the  particles 
may  be  obtained  from  these  features,  which  should  be  observable  on 
diffuse  reflection  even  against  the  background  of  the  multiply 
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scattered  field. 

d.  The  calculation  of  the  normalized  scattered  intensities 
corroborates  the  experimental  result  that,  at  least  in  the  case  of 
dielectric  particles,  the  intensity  at  a  fixed  scattering  angle 
around  45°  is  directly  proportional  to  the  total  number  concentration 
per  unit  volume  independently  of  the  size  distribution  and  size 
range.  Thus  under  certain  conditions  and  for  optically  thin  media, 
the  concentration  or  degree  of  turbidity  may  be  obtained  by  this 
method  as  veil  as  by  the  usual  one  based  on  the  extinction  of  a 


collimated  beam. 
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